Quantum spin fluctuations in the spin liquid state of Tb2Ti207 
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Neutron scattering experiments on a polycrystalline sample of the frustrated pyrochlore magnet 
Tb2Ti2 07, which does not show any magnetic order down to 50 mK, have revealed that it shows 
condensation behavior below 0.4 K from a thermally fluctuating paramagnetic state to a spin-liquid 
ground-state with quantum spin fluctuations. Energy spectra change from quasielastic scattering 
to a continuum with a double-peak structure at energies of and 0.8 K in the spin-liquid state. 
Specific heat shows an anomaly at the crossover temperature. 

PACS numbers: 75.10.Jm, 75.40.Cx, 78.70.Nx 



1. Introduction 

Magnetic systems with geometric frustration, a proto- 
type of which is antiferromagnetically coupled Ising spins 
on a triangle, have been intensively studied experimen- 
tally and theoretically for decades ^] . Spin systems 
on networks of triangles or tetrahedra, such as triangu- 
lar [3], kagome lij, spinel [HI], and pyrochlore lattices, 
play major roles in these studies. Their rich variety of 
phenomena includes zero temperature entropy in Ising 
antiferromagnets 0, Q and in ferromagnetically coupled 
spin-ice Q, multiferroics induced by non-coUinear mag- 
netic structures 0, heavy fermion behavior and un- 
conventional anomalous Hall effect ■ 

Subjects that have fascinated many investigators in re- 
cent years are classical and quantum spin-liquid states 
[llMl4| . where conventional long-range order (LRO) is 
suppressed to very low temperatures. Quantum spin- 
liquids in particular have been challenging both theo- 
retically and experimentally since the proposal of the 
resonating valence-bond state A prototype spin- 



netic LRO state at about T ~ 1.8 K within a random 
phase approximation (20| . The dynamical ground state 
is a candidate for a quantum spin-liquid, but its puzzling 
origin has been in debate [l| . 

Recently, an interesting scenario to explain this spin- 



liquid state was theoretically proposed [21|; Tb2Ti207 
is a quantum-mechanical version of the classical spin-ice 
6] , where additional spin- flip terms to the otherwise clas- 
sical Ising-spin Hamiltonian lift the macroscopic degen- 
eracy of the classical "2-in, 2-out" ground states 12 ll . 
More recently a single-site mechanism was proposed |22| 
to account for the absence of LRO, in which the CF 
ground doublet becomes two singlets owing to a conjec- 

23, 



tured tetragonal distortion 



although this inter- 
pretation is not without difficulties [25l |. 

In addition to the theoretical puzzle about the ground 
state, some experimental results of Tb2Ti207 contradict 
each other [3, 12 > HE 27 1. It was concluded that the 



liquid, or cooperative paramagnet, is the Heisenberg an- 
tiferromagnetic model on the pyrochlore lattice coupled 
by nearest-neighbor exchange interaction, which is shown 
to remain disordered at all temperatures 
ice materials, R2T2O7 (R = Dy, Ho; T 
the best-understood classical examples 



13. The spin 



Ti, Sn), are 
151, while other 



experimental candidates found recently are awaiting fur- 
ther studies 



majority of the spins is dynamic down to 50 mK in jl'^ . 
On the other hand, in 26|, |27[ it was reported that about 
50% of the spins are static at 0.4 K and that there is an 
unknown phase transition at 0.37 K. This discrepancy is 
probably caused by a certain uncontrollable parameter 
of crystalline samples. In fact, a recent study of spe- 
cific heat showed a sample dependence for single crystals 
28| . However by restricting oneself to experimental 
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13|,ll4|. 



Among magnetic pyrochlore oxides [l|, Tb2Ti207 has 
attracted much attention because it does not show any 
magnetic LRO down to 50 mK and remains dynamic 
with short range correlations [IE [l3| • Theoretical con- 
siderations of the crystal-field (CF) states of Tb^"*" and 
exchange and dipolar interactions of the system 18- 20j 



data of polycrystalline samples, results are more consis- 
tent. Muon spin relaxation (/iSR}_ 1^, neutron spin echo 
(NSE) [13, and susceptibility [29[ experiments showed 
no conventional phase transition and LRO. Only a small 
amount (~ 10 %) of spins become quasi-static at 0.1-0.3 
K 171 1221 ■ High-resolution neutron powder-diffraction 
experiments [30| showed that the crystal structure of 
Tb2Ti207 is consistent with the pyrochlore structure 
showed that it should undergo a transition into a mag- without disorder, Tb/Ti site interchange, or oxygen de- 
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ficiency. 

In this work, we hypothesize that polycrystaUine 
samples show the genuine characteristics of the spin- 
hquid state of Tb2Ti207, and reinvestigate the low- 
temperature spin fluctuations of a polycrystaUine sample 
by inelastic neutron scattering. From the NSE experi- 
ment [l7| . we expect that important spin fluctuations of 
the spin-liquid state should appear in the energy range 
E > 0.05 meV, i.e., NSE time < 0.01 ns (figure 3 of [11]), 
where no experimental data have been reported. The 
other aim of this work is to observe a certain tempera- 
ture dependence of energy spectra around ~ 0.5 K, an- 
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ticipated from the quantum spin-ice theory 2l[. In fact, 
the NSE 11 711 and our unpublished neutron scattering ex- 
periment [31I did suggest such a T dependence. We have 
found that around 0.4 K a high-temperature quasielas- 
tic spectrum becomes a continuum with a double-peak 
structure at energies of and 0.8 K, indicating that a 
crossover from the paramagnetic to spin-liquid state oc- 
curs. Specific heat shows an anomaly at the crossover 
temperature. 

2. Experiment 

PolycrystaUine samples of Tb2Ti207 were prepared by 
the standard solid-state reaction at 1350 °C from Tb407 
and Ti02 Most of the neutron-scattering measure- 
ments were performed on the triple-axis spectrometer 
NG5 at the NIST Center for Neutron Research. A sample 
with a weight of 7 g was mounted in a dilution refrigera- 
tor. The spectrometer was operated using a final neutron 
energy of Ef — 2.5 meV, providing an energy resolution 
of 0.06 meV (full width at half maximum, FWHM) at 
the elastic position. Higher-order neutrons were removed 
by cooled Be and BeO filters. A few preliminary mea- 
surements were performed on the triple-axis spectrometer 
HER at the Japan Atomic Energy Agency 3l|. Specific 
heat was measured by the heat-relaxation method on a 
physical-property measurement-system equipped with a 
•^He refrigerator. 

3. Results and discussion 

To investigate the ground state, we performed inelastic 
scattering experiments in a low energy range —0.1 < E < 
0.5 meV at low temperatures. In figures [T] and [5] we show 
constant- Q scans measured at several wave vectors in a 
range 0.6 < Q < 1.6 A-i. Since the first excited CF level 
is about 18 K or 1.7 meV [3, El], the E spectra of these 



figures are transitions among the ground doublet of the 
CF states of Tb'^"'" , being in the trigonal symmetry [3| . 
The typical energy scale of these E spectra is the order of 
0.1 meV 1 K. This value agrees well with the estimate 
of the effective Tb-Tb spin exchange-interaction based 
on Thus the observed E spectra represent 

spin fiuctuations of up and down states of the Tb^+ Ising- 
like spins [18| , interacting via the inter-site couplings and 



virtual excitations to the excited doublet 

One can see from figure[T]that the E spectra change re- 
markably between T = 1.5 and 0.1 K from a quasi-elastic 
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FIG. 1. Constant-Q scans taken at Q = 0.6, 0.8, 1.0, 1.2, 1.4, 
1.6 A"'^ in the energy range —0.1 < E < 0.5 meV measured 
at temperatures above and below 0.4 K. 
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FIG. 2. Constant-Q scans at Q = 0.8 A"^ m an energy range 
—0.1 < E < 0.5 meV measured at several temperatures down 
to 0.1 K. 



scattering centered at = to an inelastic scattering 
with an additional peaked structure at i? ~ 0.1 meV. 
Temperature dependence shown in figurc[2]indicates that 
the spectral change occurs around T = 0.4 K. We think 
that this crossover has an important implication that the 
paramagnetic states above and below 0.4 K have very dif- 
ferent characteristics. The E spectra (figure [2]) become 
a low-T limit only below 0.3 K, and we may conclude 
that Tb2Ti207 condenses into the spin-liquid ground- 
state below this temperature. It should be noted that 
in the same crossover T range a significant change was 
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observed in the NSE spectra (figure 3 of 17[). In contrast 
to this inelastic scattering, the energy-integrated diffrac- 
tion varies only slightly below 1.5 down to 0.05 K (figure 
2 of 17|), i.e., spin correlations over a single tetrahedron 
are retained at low temperatures [l6j . 

The E spectra in figure [T] above 0.2 meV at 0.1 
K exhibit additional Q-dependent structures. This Q- 
dependence seems to impose a restriction on the origin 
of the spin-liquid ground-state; it is brought about by a 



many-body effect |18l - l20l |. An interesting proposal along 
this line is the quantum spin-ice state, where a singlet 
ground state is formed predominantly from the "2-in, 2- 
out" classical spin-ice states within a single tetrahedron 
21| . The estimate of an energy band spanned by ex- 
ited states is the order 0.5 K [21|. This value approxi- 
mately agrees with the fit parameter A of the inelastic 
Lorentzian function at T < 0.3 K, which is discussed 
later. 

Quite recently another possibility of the spin-liquid 
state considering both single- and inter-site effects has 
been pointed out [22| . The authors hypothesize that the 
ground doublet in the cubic pyrochlore structure splits 
into two singlets under a small tetragonal distortion, 
conjectured to exist as high as 1.6 K [22] ■ The energy 
splitting of 0.19 meV between the two singlets, which is 
claimed to be observed at 1.6 K in a neutron inelas- 
tic spectrum in (sij . is not reproducible in the present 
data of figures [1] and [2] at 1.5 K. The energy resolution 
0.16 meV (FWHM) in [32|] seems too large to exclude 
an artifact in a resolution-convolution fitting. Thus the 
theoretical two-singlet scenario [l^] assuming the tetrag- 
onal distortion would be seriously modified to account 
for the present iJ-spectra and T-dependence, in which 
the crossover around 0.4 K could perhaps be ascribed to 
a structural transition. 

In order to parameterize the E spectra shown in figure 
[2] as a function of temperature, we carried out fits of the 
data to a scattering function 
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which is convoluted with a resolution function. The first 
and second terms of equation ([1} represent elastic and in- 
elastic scatterings, respectively. Typical results of the fit- 
ting are shown in figure|3l For data at 1.5 K we tried to fit 
the spectrum using the Lorentzian function (A = 0) for 
standard quasielastic scattering with {Aj^O, figure [3Ka)) 
and without {A = 0, figure [3Ib)) the elastic component. 
These fits show that the elastic component, 20~30% of 
the energy integrated {E < 0.5 meV) total intensity, can- 
not be neglected in this analysis. The elastic part con- 
sists of incoherent elastic scattering of Ti nuclei, multiple 
Bragg scattering, and magnetic scattering of Tb^+ spins 
17| . In the present experiment, the separation of the 
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FIG. 3. Results of fitting constant-Q scan data at T = 
1.5 K (a,b) and 0.1 K (c) shown in figure [2] to resolution 
convoluted S{Q,E) of equation ([1]). We assume (a) A = 0, 
and (b) A — and A = 0. Red dashed-lines and blue dotted- 
lines represent the resolution convolutions of the elastic and 
inelastic scattering of equation The black solid- lines are 
total fitted curves, (d) Temperature dependence of the fit 
parameters F and A. Lines are guides to the eye. 



owing to the insufficient instrumental energy-resolution 
compared to the energy scale of 0.1 meV. We could not 
obtain a reasonable T-dependence of the parameter A of 
equation ^ , and assumed the same A value for all tem- 
peratures. Resulting fits are shown in figures [3l[a) and 
(c). 

The temperature dependence of the fit parameters F 
and A are plotted in figure . These parameters show 
the crossover behavior around 0.4 K. At 0.1 K the second 
term of equation ([1} with A ~ 2F represents a quasi- 
gapped continuum with a significant spectral weight at 
E = 0. It should be noted that the very high resolution 
NSE data imply that ~ 20% of the elastic compo- 
nent in equation ([T|) contains a magnetic contribution 
below 0.3 K. Thus the magnetic spectrum, background 
subtracted equation ([T}, in the spin- liquid state proba- 
bly has a two-peak structure at E — and 0.07 meV. 
This should be confirmed more directly using another 
spectrometer with a much higher resolution ~ 10 /xeV or 
better in further work. 

We measured the specific heat Cp of the polycrystalline 
sample to check whether the crossover around 0.4 K can 
be observed in thermodynamic properties. The result is 
plotted in figure U together with previous measurements. 
Cp of the present work shows an upturn below 0.5 K, 
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FIG. 4. Temperature dependence of specific heat of poly- 
crystalline and single-crystalline samples [T^. [27l. [2^ . 



being consistent with the neutron data. We checked that 
the same upturn of Cp is observed by a polycrystalhne 
sample prepared in the same way as described in [igI . Ii7| . 
One can also see differences of Cp for crystalline sam- 
ples 3 24, 27, 2^, demonstrating the difficulty of dis- 
cussing experimental data taken on different single crys- 
tals, especially below 1 K. Control parameters of single 
crystals could be very small disorders, Tb/Ti site inter- 
change, oxygen deficiency, or local stress built-in during 
single-crystal growth carried out using image furnaces. 
We speculate that the large differences of Cp at low T 
may imply that the system is located close to a quan- 
tum critical point which is affected by these hidden ma- 
terial parameters. We hope that the mechanism of the 
spin-liquid state can be explored further in studies on 
well-characterized single crystals. 



4. Conclusion 

In summary, inelastic neutron scattering has been used 
to extend the work of [l^ 17 1 and explore the spin- 
liquid state of polycrystalhne Tb2Ti207. This system 
condenses into a quantum spin-liquid state from the para- 
magnetic state via a crossover around 0.4 K, where spe- 
cific heat shows an anomaly. The energy spectrum in 
the spin-liquid state is a continuum with a double-peak 
structure at energies of and 0.8 K. Its wave- vector de- 
pendence suggests that the spi n-liq uid state is brought 



about by many-body effects [18j, 121 
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